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Abstract 

A multi-port macromodeling method of the thin-wire 
structure system is proposed in this paper. In this method, 
the impedance and the admittance of the conductors are 
calculated using the partial element equivalent circuit (PEEC) 
method and the node voltage equation is achieved based on 
the node voltage method. In order to reduce the order of the 
model, the open-circuit suppression method is adopted. 
Making use of the vector fitting (VF) method and backward 
difference method, a time domain solution is achieved. The 
proposed method is validated by comparing its results with 
the calculated or measured ones published in other papers. 
Finally, an application about the direct lightning overvoltage 
in high-speed railway system is implemented. 
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Introduction 

The grounding grid, tower of the transmission line, 
and lightning protection system in the railway system 
can be viewed as a thin- wire system, which should be 
modeled in order to calculate the overvoltage due to 
the direct or indirect lightning striking. Many methods 
has been proposed to model the thin-wire structure 
system by using the method of moments (MOM), the 
finite-difference time-domain (FDTD) method, the 
hybrid electromagnetic model and transmission line 
approaches. Recently, the partial element equivalent 
circuit (PEEC) is adopted to analyzing the thin-wire 
structure such as grounding grid and lightning 
protection system. In some complex thin-wire 
structure such as the transmission line tower or the 
integrated grounding system in the high-speed 
railway system, the number of the conductors is very 
large which will result in a high order of the model 



and inefficient computation. In order to handle this 
problem, the open-circuit suppression method is 
adopted to eliminate the internal nodes, as a result, a 
multi-port circuit according to the external nodes is 
built. In the section II of this paper, the modified 
image method and the PEEC method is presented. 
Section III proposed a model-order reduction method 
base on the open-circuit suppression principle and the 
solution in the time domain is provided based on the 
vector fitting method combing with the backward 
difference method. In section IV, the proposed method 
is validated by comparing its simulation results with 
the ones calculated or measured in other papers. In 
section V, the integrated grounding system in the 
high-speed railway system is modeled and the direct 
lightning overvoltages of the catenary are simulated 
using the proposed method. 

Modeling Method of the Thin-wire 
Structures 

Modified Image Method 
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FIG. 1 PRINCIPLE DIAGRAM OF MODIFIED IMAGE METHOD. (A) 
CONFIGURATION OF THE MEDIUM AND THE SOURCE (B) 
THE SOURCE AND ITS IMAGE REPRESENTING THE FIELD IN 
MEDIUM 1 (C) AN IMAGE REPRESENTING THE FIELD IN 
MEDIUM 2 

A current source I x is located at a distance h from the 
interface of tow medium as shown in Fig. 1(A), in 
which the conductivity, permittivity and the 
permeability are ai , si , ui and G2 , 82 , U2, 
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respectively. A modified image current source I 2 is 
introduced and the field in the medium 1 can be 
calculated using the superposition theorem, which is 
shown in Fig. 1(B). In order to calculate the field in the 
medium 2, another modified image I 3 is introduced, 
as shown in Fig. 1(C). Both the value of the image 
current source can be evaluated according to the 
boundary condition. 

Suppose £ and s as the permittivity of the air and the 
ground, a as the conductivity of the ground, the value 
of the image current source I 2 and / 3 can be evaluated 
using the following equations. 

Type I: Current source and the observation point both 
in the soil, the field can be evaluated as a sum of the 
field of the current source I x and its image I 2 : 



j JCOS + CJ- JCO£ Q . 
2 j(D£ + <J + jC0£ Q 1 



(1) 



Type II: Current source and the observation point both 
in the air, the field can be evaluated as a sum of the 
field of the current source I x and its image I 2 : 



j = jcos {) - jcos-a . 
2 jco£ + <j + ja>£ 1 



(2) 



Type III: Current source in the soil and the observation 
point in the air, the field can be evaluated as field due 



to the modified current source I 3 : 



2jO)£ Q 



jC0£ + <J+ jC0£ Q 



(3) 



Type IV: Current source in the air and the observation 
point in the soil, the field can be evaluated as field due 
to the modified current source A : 



3 ~ • , , • 1 

JC0£ + (J + JCO£ {) 



(4) 



PEEC Method 



PEEC method is a full wave method for solving the 
conductor system based on electric-field integral 
equation, and a lumped circuit model can be obtained. 
Based on the principle of this method, the thin-wire 
structure can be viewed as a combination of short 
filament whose current density and charge density 
distribute along the conductor. The total tangential 
electric field shall be as (5). 



s • E* (? ) = s • [i? (? ) + E 5 (? )] = s 



(5) 



where, J is a unit tangential vector along the wire, 
E l (?) is the incident electric field, E s (?) is the scattered 
electric field which can be represented as (6). 



E s (r) = ja>A(r)-V<j)(r) 



(6) 



where, A(?) is the magnetic potential vector and ^(?) 
is the electric scalar potential, whose expression 
equation is 

A(r) = -^- [?./(?).-* 



47ri 
1 



R 



>ds* 



kr) = -—\p l ir)' e --'ds^ 
4x£ J R 



(7) 
(8) 



where, ds ' is a element of the wire, T = yJjcoju(<j + jcc>£) 
is the propagation constant of the considered medium, 
the distance from the source and the observation point 
isR = |?-?'| . 

Substituting (6), (7) and (8) into(5), (9) is obtained 



s-E l (r)-s 



iCr) >~»-\s-mr) e ' 



--jco- 



rR ds'-*l = (9) 



R ds 

where a l is the conductivity of the conductor. 

Considering a segment connecting with node k and I, 
and Integrating (9) along the segment from, (10) is 
obtained. 



s ' =l P - TR 
k=h-k=-j«>-T-_ \ \s-s'I(r') e 



4tt ? , 
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(10) 



According to (10), the first term can be defined as the 
impedance of the along the conductor segment, which 
can be expressed as 



^Lmn 



An 



e FR cos cp 



R 



-dldl 



(11) 



Based on the conservation of charge, the charge 
density t, can be expressed as 



(12) 



where, I Ti is the leakage current and the L k is the 
length of the fcth segment. As a result, the average 
potential can be calculated using the following 
equations. 



4^(cr+ jco£)L k / J c L R 



\\f* e —dldl (13) 
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Let 

M 

(4=E Z ™4 (i4) 

z=l 

The horizontal impedance Z m can be defined as 

Z m = . — — \\—dhdl k (15) 

As we can see, the form of Z m is the same with Z Lmn , 
which can be calculated using Neman integration or 
numerical method. 

Multi-port Macromodel of the Complex Thin- 
wire Structures 

Model-order Reduction Method based on the Open- 
circuit Suppression Method 



A thin-wire conductor system can be viewed as a 
number of cells composed of straight filaments, whose 
PEEC circuit is shown as Fig. 2. 




FIG. 2 THIN-WIRE STRUCTURE AND ITS EQUIVALENT CIRCUIT 

According to a certain filament k connecting between 
the node i and j, whose electric potentials <j> k can be 
expressed by (16) because the length of the filament is 
short enough. 



h=^ L Y JL ( 16 ) 

where, U ni and U nj is the node voltage of node i and j, 

respectively. Considering all of the filaments in the 
conductor system, the relation between the filament 
electric potential and the node voltage can be written 

as 

(j) = 0.5 x abs( A T ) U n = P U M (17) 

where, A is the incidence matrix in the node voltage 
method, the function abs(-) denotes to a calculation of 
the element's absolute value in the corrsponding 
matrix. 

Based on the node voltage method, the node 



admittance matrix considering only the axial direction 
impedance can be expressed as (18). 

Y^=A.V-A r (18) 

where Z b axial direction impedance. 

Meanwhile, the leakage current can be described by 
the leakage admittance matrix and the relation 
equation between the filament electric potential and 
the leakage current is 

i»=Y»^ (19) 

where, is the leakage admittance matrix, which can 
be obtained by = Z^ 1 . 

Assuming the leakage current flowing from the two 
nodes of each filament which is equal to each other, 
the relation equation of the leakage current and the 
node current can be written as (20) 

l n =0.5xabs(A)-i b =Q l b (20) 

Based on the Kirchhoff's circuit laws, the equations of 
the circuit can be obtained, as shown in (21). 

-1^=1,-1, (21) 

where, l s denotes to the incident current vector of the 
current source. Substituting (18), (19) and (20) into (21), 
the equations corresponding to the node voltage 
method is 

(A.Z^A r +Q.Y,Q r ).U„=i 5 (22) 

Let Y n = A Z, 1 A r +Q Y,Q r , and Y n is the node 
admittance matrix of the whole network. 

For a thin-wire system with a large number of 
conductors, the order of the node admittance matrix 
would be very large which would result in an 
inefficient computation. In this section, the open- 
circuit suppression principle is used to reduce the 
order of the model. 

The nodes in the conductor system can be divided into 
two types according to the different location of the 
nodes, one kind is the node connecting to the external 
circuit and another is internal nodes. Because the 
internal nodes would not connect to the external 
circuit, they can be eliminated to reduce the order of 
the model. 

Let the voltage and current vector of the internal 
nodes are U 2 and i 2 , the voltage and current vector of 
the internal nodes are U 1 and I 1 . Rearrange (22) as 
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(23) 



Considering the incident current of the internal nodes 



: , so the second equation of (23) can be written 



as 



Y 21 U„ 1+ Y 22 U„ 2 



(24) 



and the node voltage of the internal nodes can be 
expressed by the node voltage of the external nodes, as 
shown in (25). 



U =-Y 1 Y U 



(25) 



Substituting (25) into the first equation of (23), the 
equations of the external nodes can be obtained. 



I, 



:(Y -Y Y _1 Y )t\ 
\ A ll A 12 A 22 A 21 J^nl 



(26) 



From equation (26), a new node admittance matrix is 
obtained, as shown in (27). Because the number of the 
external nodes is much less comparing with the 
number of the total nodes, the order can be reduced 
effectively. 



Y' = Y„ 



Y Y Y 

x 12 x 22 x 21 



(27) 



Time Domain Solution of the Macr onto del 

For the purpose of simplicity, a simple grounding 
system is used to represent the mathematical 
derivation of the solution of the macromodel, as Fig. 3. 
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5 -5 

FIG. 3 SIMPLIFIED GROUNDING SYSTEM 

Step 1: Based on the method of previous section, the 
node admittance matrix can be obtained and it can be 
fitted using the vector fitting method, as shown in (28). 
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(28) 



Step 2: Using the Laplace inverse transformation, the 
time domain equation can be obtained, shown in (29) 



a { 
a t 



d l u x 
df 
d l u 2 
. df 



^11 i + ^11 " a i \l i + ^12 " a i 

^21 i + d n ' a i b 22 . + d n • a. 



d l i x 
df 
d l i 2 
~df 
(29) 



Step 3: Making use of the backward differentiation 
formula to approximates the time derivatives, the 



I 

i=0 

=-X 

where, 



a t 
a ; 

NF 



~e^,\ ; c n 



21,i ^22,i 

a r B a (i,:) -c Ui -c l2 . 
a r B a (i,:) -c 2U -c 22 . 



(30) 



B = 



co NF (NF) a> NFl (NF) co NF _ 2 (NF) 
cd nfa (NF-1) m NF _ 2 (NF-\) 















a\(NF) 



, CL)j(p) = (-l) J .J , which is a cumulative sum 
function, x his =[u* ls u 2 his if i 2 his J , in which 

u His = [ u n~NF + l ... ^ ^ ^his = ^n-NF+l ... R n j ^ 
{ k lS= ^n-NF + X ... .„j f ^s^n-NF+X ... ^ ^ ^ 

\ n+1 = [u x n+1 u 2 +1 i x n+l i 2 +l ~^ . Combined with the 

boundary conditions of the model, a transient analysis 
can be accomplished. 

Validation of the Method in Thin-wire 
Structure System Transient Analysis 

Transient Analysis of a Vertical Conductor above the 
Ground 
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FIG. 4 CALCULATIONAL MODEL OF A VERTICAL 
CONDUCTOR INJECTED AT THE TOP BY A STEP CURRENT 
SOURCE. (A) CALCULATION MODEL (B) INJECTED CURRENT 
SOURCE 

The calculation model consists of a 60m high vertical 
conductor connected to a grounding resistance of 10Q 
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at its bottom and to a surge impedance of 1500Q at its 
top, which is shown in Fig. 4(A). A step current was 
assumed to be injected at the top of the conductor. 
Three different front times were considered, namely tf 
=0.2 us, f/=0.5us and £/=lus, as shown in Fig. 4(B). 

Transient voltages at the top of the vertical conductor 
according to different front times are calculated using 
the macromodel proposed in this paper, as shown in 
the left side of Fig. 5, and the corresponding results 
published in Alberto D. Conti's paper are shown in the 
right side. It is shown that the all of the results match 
well with each other, which demonstrates the accuracy 
of this method for dealing with the transient analysis 
of the vertical conductor. 



- HEM 

- ATP (ore section) 

_ AtP(inree sections) 




FIG. 5 VOLTAGES AT THE TOP OF THE VERTICAL 
CONDUCTOR. (A) RESULTS ACCORDING TO £/=0.2|as (B) 
RESULTS ACCORDING TO tf=0.5[is (C) RESULTS ACCORDING 
TO tf=l\is. 

Transient Analysis of a Vertical Conductor System 
above the Ground 

A vertical conductor system of four parallel vertical 
conductors with 16.5 mm radium is shown in Fig. 6. 
This is the same configuration carried out by Hara. 
The vertical conductors are excited by a pulse 
generator (PG) through a current lead wire, and the 
top voltage is defined as the voltage between the 
conductor top and a measuring wire. In the simulation, 
the resistivity of the earth was set to 1.69xl0 8 Q-m, 
PG was modeled by a current source with 5kQ 
internal resistance, whose waveform is shown in Fig. 7. 




FIG. 6 CONFIGURATION OF A VERTICAL CONDUCTOR 
SYSTEM EVALUATED BY HARA 
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FIG. 7 WAVEFORM OF THE CURRENT SOURCE 




10 20 30 40 
Time(ns) 



FIG. 8 VOLTAGE AND CURRENT CALCULATED AND 
MEASURED AT THE TOP OF THE VERTICAL CONDUCTOR 
SYSTEM. (A) CALCULATED RESULTS (B) MEASURED BY HARA 

The waveforms of the tower-top voltage and current 
are shown in FIG. 8 (A). The results measured by Hara 
are given in Fig. (B). Both the results agree well with 
measured ones which demonstrate the validity of the 
macromodel proposed in this paper. 

Transient Analysis of a Grounding Grid 

A 60 x 60m 2 grounding grid with 6 by 6 10 meter square 
meshes, 10m ground rods in the corners constructed of 
copper conductors with diameter 1.4cm, is buried at 
0.6m depth, as shown in FIG. 9. The soil's resistivity 
and relative permittivity is 100 Q-m and 36, 
respectively. A 1/20 us lkA crest current surge is 
injected at the middle point of the grid. 
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FIG. 9 GROUNDING GRIDS IN THE SIMULATION 
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Transient voltages at the feeding point A, middle point 
B and corner point C are presented in Fig. 10(A). 
Results calculated by Grcev Leonid D. are illustrated 
in Fig. 10(B). The good agreement between the 
compared results has been achieved. 




- Point A 

■ Point B 

■ Point C 



10 
Tirme(|js) 




FIG. 10 COMPUTED TRANSIENT POTENTIAL AT THREE 
POINTS ON THE GROUNDING GRID (A) POTENTIAL OF A, B 
AND C CALCULATED BY PEEC (B) RESULTS CALCULATED 
GRCEV LEONID D. 



grounding sticks at the corner the pile cap are adopted 
as the grounding rod and their length are 5~30m in 
general. The cross section area of the reinforcing bars 
should not less than 200mm 2 . 



Axial direction grounding steel bar 



Connection point between the 
pier and underground 
grounding systei 




Connection point between the 
box girder and the pier 



Applications 

Because the viaduct structure are widely adopted in 
the China's high-speed railway system, the increase of 
the height of the catenary's overhead lines will lead to 
a frequent lightning striking. Nowadays, there is no 
ground wire equipped above the power supply lines 
whose geometry is shown in Fig. 11, which may result 
in a serious lightning accident. In order to represent 
the effect of the ground wire and evaluate of the 
withstand level of the catenary, the integrated 
grounding system is modeled and the lightning 
over voltage of the catenary is calculated. 




FIG. 11 GEOMETRY OF CATENARY IN THE HIGH-SPEED 
RAILWAY SYSTEM 

Structure of the Integrated Grounding System of the 
Viaduct in High-speed Railway 

The structure of the integrated grounding system of 
the viaduct is shown in Fig. 12, in which there are 6 
grounding reinforcing bars are located on the surface 
of the box girder and 2 vertical grounding reinforcing 
bars are use to connect the box girder and the pier, 
meanwhile, the structural reinforcing bars can also be 
used as the grounding conductors in principle. The 



FIG. 12 STRUCTURE OF THE INTEGRATED GROUNDING 
SYSTEM OF THE VIADUCT 

Macromodel of the Integrated Grounding System of 
the Viaduct in High-speed Railway 

Because of the complex structure of the integrated 
grounding system, it is divided into three sections in 
this paper, they are the grounding system of the box 
girder, the grounding system of the pier above the 
ground and the grounding system underground, as 
shown in Fig. 13. 




5 5 



FIG. 13 MODELING OF THE INTEGRATED GROUNDING 
SYSTEM OF VIADUCT 

Using the macromodeling method proposed in this 
paper, the multi-port macromodel of each section can 
be obtained and the whole model can be built 
according to the connection point of the integrated 
grounding system, which is shown in Fig. 14. 
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Multi-port model of box 
girder grounding system 




Multi-port model of 
pier grounding 
system 



FIG. 14 EQUIVALENT NETWORK OF THE INTEGRATED 
GROUNDING SYSTEM OF VIADUCT 

The overhead lines of the catenary can be modeled as 
a multi-conductor transmission line system, combined 
with the multi-port macromodel of the integrated 
grounding system, the lightning overvoltage is 
calculated due to a lkA 2.6/50 [is lightning current 
directly striking at the ground wire. The voltages of 
the power supply line's insulator are shown in Fig. 15, 
in which the resistivity of the ground is 100 Q-m, and 
the withstand level of the catenary are listed in Table 1 
with different ground resistivity. 
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FIG. 15 POTENTIALS OF THE AF LINE AND T LINE AND 
VOLTAGES OF THE INSULATOR OF AF LINE AND T LINE 
(HEIGHT OF THE VIADUT IS 8M, DEPTH OF THE GROUNDING 

STICK IS 8M AND EARTH RESISTIVITY ( p = 100Q • m ). (A) 
POTENTIALS OF THE AF LINE AND T LINE (B) VOLTAGES OF 
THE INSULATOR OF AF LINE AND T LINE 

Table 1 Lightning withstand level of AF line and T line due to 

LIGHTNING STRIKING ON GROUND WIRE 



Resistivity of 
the ground 
(n-m) 


100 


500 


1000 


Catenary's line 


AF 


T 


AF 


T 


AF 


T 


Withstand level 
(kA) 


58.79 


58.13 


58.50 


57.83 


58.49 


57.81 



The lightning withstand levels of AF line and T line 
show little difference with different ground resistivity, 
that is because the large number of the grounding 
conductors resulting in a good grounding effect. 

In order to represent the effect of the ground wire, the 
overvoltage due to a direct lightning striking to AF 
line and T line are also calculated and the withstand 



levels are shown in Table 2. It can be concluded that 
the equipment of the ground wire can reduce the 
overvoltage effectively, which will improve the 
reliability of the power supply in the high-speed 
railway system. 

Table 2 Lightning withstand level of AF line and T line due to 

LIGHTNING STRIKING ON AF LINE AND T LINE 



Catenary's line 


AF 


T 


Withstand level (kA) 


1.79 


3.49 



Conclusions 

A time domain multi-port macromodeling method of 
the complex thin-wire system is proposed in this paper. 
Based on the method, a circuit with low order can be 
achieved and a time domain numerical solution 
method is provided which can ensure the convergence 
of the model. The accuracy of the proposed method has 
been validated by comparing the results with the 
simulation or measurement results published in other 
papers. 
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